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Working in the linear sigma model with quarks, we compute the finite-temperature effective 
potential in the presence of a weak magnetic field, including the contribution of the pion ring 
diagrams and considering the sigma as a classical field. In the approximation where the pion self- 
energy is computed perturbatively, we show that there is a region of the parameter space where 
the effect of the ring diagrams is to preclude the phase transition from happening. Inclusion of the 
magnetic field has small effects that however become more important as the system evolves to the 
lowest temperatures allowed in the analysis. 
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I. INTRODUCTION 

In recent years it has been possible to produce and 
study hadronic matter at high densities and temperatures 
by means of collisions of heavy nuclei at high energies [l[ . 
There are convincing signals that reveal the production 
of deconfined matter where the degrees of freedom in- 
volved are the quarks and gluons of QCD, forming the 
so called quark-gluon plasma (QGP). A common feature 
of these signals is their strengthening as the centrality of 
the collision increases. 

At the same time it has been realized that a host of new 
phenomena can also happen for not so central collisions. 
Among these, it has been pointed out that for peripheral 
collisions, a magnetic field of a non-negligible strength is 
generated Q. The origin of this field is two- fold: On one 
hand, in non-central collisions, there is a local imbalance 
in the momentum carried by the colliding nucleons in the 
target and projectile that generates a non-vanishing local 
angular momentum [H, Q which in turn produces a mag- 
netic field, given the net positive charge present in the 
collision. On the other hand, the spectator nucleons can 
be thought of as currents of net positive charge moving 
in opposite, off center, directions which in turn produce 
a magnetic field that adds up in the interaction region. 

An interesting question that emerges from this scenario 
is whether a magnetic field can influence the phase tran- 
sitions that may occur during the reaction, in particular 
the chiral phase transition. There are known examples 
where magnetic fields are able to change the nature of 
a phase transitions. Most notably is the Mcissncr ef- 
fect where the phase transition of type I superconduc- 
tors changes from second to first order when an external 
magnetic field is applied. Magnetic catalysis is another 
phenomenon whereby a magnetic field is able to dynami- 
cally generate masses in QED, regardless of the strength 
of the field Q. More recently, it has been shown that 



in the presence of primordial magnetic fields, the elec- 
troweak phase transition, that took place in the early 
universe for temperatures of order 100 GeV, gets also 
strengthened Q . In a similar connection, the dynamical 
generation of anomalous magnetic moment of the elec- 
tron has also be unveiled in Ref. Q- 

Calculations of the intensity of the field produced in 
this kind of collisions show that for very early proper 
times after the reaction (r < 0.1 fm) the field reaches 
values eB ~ 6(m™ c ) 2 , where to™ c is the vacuum pion 
mass, even for mid-peripheral collisions. The intensity 
decreases with the proper time r as eB oc 1 /t 3 in such a 
way that for t ~ 1 fm, namely, for times when the stan- 
dard picture of a heavy-ion reaction places the existence 
of the equilibrated QGP, eB < 0.1 (m™ c ) 2 0, that is, 
already two orders of magnitude smaller than at the very 
early stages of the collision. 

In a recent work, the chiral phase transition in rel- 
ativistic heavy-ion collisions has been examined in the 
presence of strong magnetic fields using the linear sigma 
model @. Working with the hierarchy of scales where 
(ra™ c ) 2 < T 2 < eB, with T being the temperature 
around the phase transition, the authors conclude that 
the effect is to turn a crossover into a weak first order 
transition, a result that might be relevant even for the 
physics of the primordial QCD transition. Nevertheless, 
as mentioned above, a more realistic scenario in heavy- 
ion collisions should be to consider that for the times 
when the initial chromoelectric fields decohere in the af- 
termath of the collision and give rise to partons -which 
in turn are the appropriate degrees of freedom to de- 
scribe the chiral phase transition the magnetic field in 
the interaction region might not be that strong. At these 
times the hierarchy of scales is such that the magnetic 
field is the smallest of all and the deconfinement/chiral 
phase transition temperature is the largest one, while 
the pion vacuum mass, occupies an intermediate place. 
Furthermore, the analysis of Ref. @ neglects the contri- 
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bution from the so called ring diagrams which arc known 
to be important at high temperatures to account for the 
infrared properties of the plasma 0, licj . Similar con- 
siderations, in QED for instance, have already been re- 
ported ll| . In the context of the QGP and hadron mat- 
ter, the effects of magnetic fields have also been recently 
looked at. Such studies include their influence on making 
evident possible topological-charge transitions in heavy- 
ion collisions [1, [Hi, on the confincment/deconfinement 
phase transition -in the Abelian approximation of the 
chromomagnetic field- [l3| and on the hadron struc- 
ture Q. 



In this work we undertake the calculation of the effec- 
tive potential at finite temperature, in the presence of 
weak magnetic fields. We use the linear sigma model as 
the working tool to describe the chiral phase transition in 
relativistic heavy-ion collisions. Wc work explicitly with 
the hierarchy of energy scales where eB <C to 2 <C T 2 , 
with m being a generic mass appearing in the calcula- 
tion and considering that the interaction region is subject 
to an external magnetic field directed along the positive 
z-axis. We work up to the contributions of the ring di- 
agrams since for the aforementioned hierarchy of scales, 
the effects of the magnetic fields appear only at this level. 
Computation of the ring diagrams calls for the calcu- 
lation of the pion finite temperature self-energy in the 
presence of the field. We approximate this self-energy by 
its perturbative value although, for the large couplings 
involved, the calculation most likely overshoots the true 
result. Wc parameterize the lack of an accurate non- 
perturbative treatment introducing a parameter to con- 
trol the strength of the self-energy and analyze the con- 
sequences in this parameter space. To incorporate the 
weak magnetic field we use the Schwinger proper time 
method to write the charged particle propagators in a 
weak field expansion. 



The work is organized as follows: After a brief sum- 
mary in Sec. |TT] of the linear sigma model Lagrangian, 
we perform the calculation of the finite temperature ef- 
fective potential in Sec. IIII1 considering that the sigma 
field is classical and going up to the contribution of the 
ring diagrams in the presence of a weak magnetic field. 
In Sec. IIVI we present the numerical results of the work 
and show that in a certain region of the parameter space, 
even in the absence of magnetic fields, the ring diagrams 
preclude the development of the phase transition. Inclu- 
sion of the weak magnetic field has little impact on the 
phase transition. Finally in Sec. [V] we discuss our results 
and present the conclusions. 



II. THE LINEAR SIGMA MODEL WITH 
QUARKS 



The Lagrangian for the linear sigma model is given by 
£ = \(d^) 2 + l(d^) 2 + ^(a 2 +n 2 ) 
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igi>T'Y5i> ■ 7T - griper, 



(1) 



where ip is a SU(2) isospin doublet of massless quarks, 
7r = (-7Ti, 7T2, ^3) is an isospin triplet representing the pi- 
ons and a is an isospin singlet. 

When the mass parameter fi 2 is positive, the La- 
grangian admits a broken symmetry vacuum solution 
given by the minimum of the classical potential 
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(2) 



Choosing this minimum along the a direction, the vac- 
uum expectation values for the sigma and pion fields are 
given by 



v = fi/\/X 
0. 



(3) 



vq is also called the classical vacuum which minimizes 
the action for uniform field configurations. To study the 
quantum properties of the system, we define the shifted 
field a' by 

o- = v + o-', (4) 

where v is taken as a variable. When v = vq, cr' represents 
the field configuration around the classical vacuum. 

As a function of the shifted field, after symmetry 
breaking, the Lagrangian of Eq. (fT]) becomes a theory de- 
scribing a a' field, three pion fields and a quark-doublet 
field with masses given by 



m 2 ,(u) 
m\(v) 
m q (v) 



3Xv 2 

Xv 2 - 

gv. 



/' 
,,2 



(5) 



Armed with the fundamentals of the linear sigma model, 
wc proceed to compute the effective potential. 



III. THE EFFECTIVE POTENTIAL 



A. Tree and one-loop 

The tree level potential is given by 
V 2 2 , A 



y(tree) 



-V + -V 1 . 



(6) 
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Wc consider that the a' field is very heavy and thus treat 
it only classically. To one-loop order, the contribution to 
the effective potential in the imaginary-time formulation 
of thermal field theory is given, for bosons, by 

= ^E/t^M^- 1 ) 1/2 ) (7) 



V, 



(i) 



whereas for fcrmions, by 



(2tt) 3 



(8) 



where Sbj are the degeneracy factors accounting for 
the internal degrees of freedom for bosons (isospin) and 
fcrmions (isospin and color), respectively, n is the index 
for the Matsubara frequency and D and S represent the 
boson and fermion Matsubara propagators. For charged 
particles, these propagators should include the effect of 
the external magnetic field. Nevertheless, it has been 
shown Q that for the hierarchy of energy scales consid- 
ered, the terms containing the effects of the magnetic 
field are subdominant. 

Equations (JT]) and (jSJ) contain both a vacuum and a 
finite temperature piece. The vacuum piece exhibits the 
usual ultraviolet divergence that needs to be removed. 
In the present context, this means that, up to additive 
constants, the effective potential should contain only fi- 
nite independent terms. This is accomplished, for in- 
stance, by introducing counter-terms to absorb the in- 
finities. The usual physical conditions implemented to fix 
the counter-terms require that the position of the mini- 
mum of the effective potential, as well as the mass of the 
a' field maintain their classical values [ll|. However, for 
theories with massless modes -such as the pions in the 
present case, whose mass vanishes at v = Vo~ this proce- 
dure breaks down. The problem is that in order to get 
the mass of the a' field from the effective potential one 
requires computing the second derivative, since this mass 
is the inertia along the c'-axis. This derivative turns out 
to not be defined at such value. Thus, the second condi- 
tion needs to be replaced by another appropriate one in 
a manner that we proceed to explain. 

First, let us introduce the general expression for the 
one-loop renormalized effective potential, where we add 
the counterterms to absorb the ^-dependent infinities 
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o(A) - Sn 2 J 2 
1 v 



b(A) + 6X\ 4 OT , 

J i J 4 + 3/(m 7r ,A) 

247(m/,A). 



(9) 



The last two terms account for the pion and fermion con- 
tributions to the vacuum effective potential, respectively, 
which from Eqs. Q and (JSj) involve the ultraviolet cutoff 
(A) dependent function I(m, A) defined as 



1 /-A 

J(m,A) = / dkkWk 2 +: 
27r Jo 



(10) 



The counter-terms in Eq. ([9]), a(A) and 6(A), are intro- 
duced to take care of the w-dependent infinities, whereas 
the counter-terms S/j, 2 and <5A, account for finite terms 
that might shift the coefficients of the v 2 and i> 4 terms, 
respectively. 

After some straightforward algebra where the in- 
dependent infinities arc absorbed, the one-loop renormal- 
ized effective potential becomes 
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(11) 



where, hereafter, when not explicitly indicated, the pion 
and quark masses are the w-dependent ones, given in 
Eqs. ©. 

To fix one of the remaining counter-terms, either <5/z 2 
or 5X, we impose the condition that the minimum of the 
renormalized effective potential remains at its classical 
value, namely, 

1 d 



2v dv ren 



0, 



which implies that 

.4 



5X = 6 



9 
4tt 2 



1 + ln 



4A 



A A 2 



(12) 



(13) 



Inserting Eq. (fTB"]) into Eq. (fTT|) . we get 



vS}} = i V 



- 6 




^ln(^ 
64tt 2 V 4 



(14) 

To fix the second counter-term, <5/^ 2 , notice that the ar- 
gument of the last logarithmic function is dimcnsionfull, 
though the whole term is well defined as m r goes to zero. 
However, if we choose 



Sfi 2 



we obtain 



~ 3 16^ ln 



A 



2^ 



64tt 2 ' 
A 2 



(15) 
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32tt 5 
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16tt 2 



In 



m 2 Jv ) 



64tt 2 



hi 



(16) 
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which no longer contains dimension-full arguments of the 
logarithmic functions. The choice of Eq. (fT5|) although 
seemingly arbitrary, has the advantage of producing a 
well defined renormalized effective potential which pre- 
serves the properties of the tree level one, namely, that 
the pions are massless at vq which in turn keeps being 
the minimum of the potential. Recall that the effective 
potential is not in itself an observable; only physical prop- 
erties extracted from it, such as the position of the min- 
imum and the critical temperature are. The choice that 
produces Eq. ([TB]) is also extensively used in Standard 
Model calculations @. 

We now proceed to include the finite temperature con- 
tribution at one-loop. The finite temperature pieces of 
Eqs. ([7]) and ||HJ) are given by 
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(17) 



Adding the renormalized effective potential to the finite 
temperature contributions, we get the full one-loop finite 
temperature effective potential 
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A few words about the properties of Eq. (|18|) are in order: 
First, notice that the dependence on the pion mass in the 
argument of the logarithmic functions has canceled upon 
addition of the vacuum and finite-temperature pieces of 
the renormalized effective potential. This is an impor- 
tant property for otherwise, this function can develop an 
imaginary part when the pion mass is negative, namely 
for v < i>o . Secondly, notice the appearance of a cubic 
pion mass term. This is also a dangerous term since it 
gives rise to an imaginary piece when the pion mass is 
negative. However, as we will show, this term is exactly 
canceled when considering the contribution from the ring 
diagrams. 

In order to closely examine the behavior of the renor- 
malized finite-temperature effective potential with the 
temperature, let us re-express Eq. (|18[) expanding it in 



powers of v 
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Let us ignore, for the time being, the term propor- 
tional to m%. Notice that the critical temperature T c for 
the phase transition is determined by the curvature of the 
effective potential at v = 0. When the curvature changes 
sign the phase transition starts. At high temperature the 
effective potential at v — is convex and the minimum 
of the potential happens for v = 0, that is the symmetric 
phase. However, when the effective potential becomes 
concave at v = 0, the minimum of the effective potential 
is located at a finite value of v, that is the broken symme- 
try phase. The above takes place provided the coefficient 
of the term v 4 is positive for otherwise the effective po- 
tential for large v becomes concave and thus unstable. 
This can happen for very high temperatures larger than 
a given temperature T max . Therefore, the condition for 
the analysis to be valid is that the critical temperature is 
smaller than this last temperature, namely, T c < T max . 

To test whether the analysis is consistent, we proceed 
to compute these temperatures. For this purpose, we use 
the standard values for the parameters 

A = 20 

fx = 380 MeV 

g = 3.3, (20) 



(18) which are determined from requiring that the vacuum a 
mass is, m v a ac ~ 600 MeV, the constituent quark mass 



m vac ^ 3qq MeV and by using a value for the pion vac- 
uum decay constant f n = 93 MeV. The critical tempera- 
ture is computed from solving for the temperature where 
the second derivative of the effective potential with re- 
spect to v vanishes, the maximum temperature is com- 
puted by finding the temperature for which the coefficient 
of the quartic term vanishes. This gives 



T c = 147.5 MeV 
T max = 7, 196.5 MeV, 



(21) 



which confirms that T c < T max . The phase transition be- 
comes a smooth crossover. Figure [1] shows the potential 
Vren , dropping out ^-independent terms, scaled by 
T 4 , for different values of the temperature and ignoring 
the cubic term in the pion mass. Notice how this poten- 
tial flattens out continuously as the temperature lowers 
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FIG. 1: Wen scaled by T c 4 and ignoring ^-independent 
terms as well as the cubic term in the pion mass, for several 
values of the temperature. 



down up to the temperature where the phase transition 
takes place, after which, the potential develops a mini- 
mum at a finite value of v. 



B. Ring diagrams 

We now proceed to include in the analysis the contri- 
bution from the ring diagrams. As is well known, for 
theories with massless modes, pcrturbativc calculations 
can lead to the appearance of infrared divergences which 
signal the need of a resummation scheme. The leading 
divergences can be summed up to render an infrared safe 
quantity and the diagrams corresponding to this diver- 
gences are known as the ring diagrams. These are de- 
picted in Fig. [2] and the explicit expression for the re- 
summed series is given by 

y(ring) _ _ T f d 3 k 

2 (2tt)3 

oo 1 

x ]T^{2[n*A B (fcr + [n AW}, 



N=l 



d 3 k 



2 ^ J (2tt) 3 

x |21n[l + FT A B (fc)] + ln[l + n°A°(fc)]}, 

(22) 

where II* and II are the charged and neutral pion self- 
energies in the presence of the magnetic field and A B and 
A their corresponding propagators. The factor 2 takes 
into account that there are two charged pions. 

At this point we should mention that, as the couplings 
(expansion parameters) in the expression for the pion 
self-energy are too large for a pcrturbativc calculation to 
be strictly valid, it is likely that the one loop expansion 
overshoots the exact result. For the purposes of this work 



FIG. 2: Schematic representation of the resummation of the 
ring diagrams. 



where we explore whether a weak magnetic field may have 
some influence on the dynamics of the phase transition, 
we will content ourselves with a qualitative description 
of the problem. For this, we parameterize our ignorance 
on the self-energies computing them perturbatively up 
to one-loop order but multiplying these times a constant 
< c < 1 that we will vary to explore the parameter 
space. The correct line of action is to calculate the pion 
self energy non perturbatively. The continuum approach 
in this connection would be to solve its Schwinger-Dyson 
equation consistently in the presence of a heat bath and 
uniform external magnetic field, and then substitute the 
result into Eqs. (|2"2"|) . However, in the present context, we 
replace this procedure by introducing the free parameter 
c as mentioned before. 

Equation (j2"2")) contains both vacuum as well as In- 
dependent infinities. These last can be canceled by a 
one-loop vacuum counterterm that renormalizcs the pion 
mass. The procedure is best carried out by explicitly 
separating the two-loop contribution to Eq. (f2"2"| which 
results in the expression [l6| 



y(ring) _ ±_ 

2 ^ 



d 3 k 



{2ir) 3 

x {2(ln[l + n*A B (fc)] -FTA s (fc)) 

+ (ln[l + n°A°(/c)] -n°A°(fc))} 

+ 2V* {2) +y 0(2) , (23) 

where V*^ and V ^ are the contributions to the two- 
loop effective potential for a charged and a neutral pion. 
Once again, the factor 2 takes into account that there are 
two charged pions. 

For the charged pion contribution to the one-loop pion 
self-energy, we use the charged scalar propagator which 
in the weak field limit is given by 
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,(24) 



k 2 - 



where u> n = 2nirT, with n an integer, E 
with k 2 ^ — k\ + k 2 and the u-dependent pion mass. 
For the neutral pion contribution to one the loop pion 
self-energy, we use the same propagator in Eq. (|2"4"]) set- 
ting eB = 0. The quark contribution to the pion self- 
energies should in principle consider that quarks are sub- 
ject to interact with the external magnetic field. Never- 
theless, for the hierarchy of energy scales we work with, 
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FIG. 3: Effective potential scaled by Ti(0.25) 4 , ignoring v- 
independent terms, computed for A = 20, fj, = 380 MeV and 
g = 3.3 in the absence of a magnetic field, for several temper- 
atures. 



FIG. 4: Effective potential scaled by 7i(0.38) 4 , ignoring v- 
independent terms, computed for A = 20, [i = 380 MeV and 
g = 3.3 in the absence of a magnetic field, for serveral tem- 
peratures. 



it is easy to see that the independent part of this self- 
energy is subdominant and thus we just consider the B- 
independent quark propagator. 

The dominant contribution in Eq. (|23[) comes from the 
mode with n = 0. Fcrmions do not contribute to the 
ring diagrams as their mode with n = does not vanish. 
It is easy to check that, after mass renormalization, the 
computation of Eq. (|23[) (for the mode with n = 0) re- 
duces to considering only the T-dependent terms. As it 
is outlined in the appendix, carrying out an expansion of 
the argument of the logarithms and keeping only terms 
up to O (eB) 2 , the dominant contribution from the ring 
diagrams, according to the hierarchy of energy scales we 
are considering, can be written as 
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(25) 

where in the fourth term on the right-hand side of 
Eq. (|25|) we have kept only the leading contribution in 
(eB). To cope with the lack of a non-perturbative cal- 
culation for the self-energies, from now on we set IT , II* 
and 111 as given by 
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(26) 



with < c < 1 and where, as sketched in the appendix, 
we make use of the same cancellations that effectively 
substitute m T — > m T , with 



(27) 



as in the calculation of the ring diagrams. There are sev- 
eral properties of Eqs. (|25p and (|26p worth mentioning: 
First, notice that, as promised, the cubic pion mass term 
in Eq. ([25| exactly cancels the one appearing in Eq. (fl9|) . 
Second, as outlined in the appendix, the leading contri- 
bution in (eB), in the first and second of Eqs. (|2"6"|) need 
only to consider II i as a correction to the pion mass. 

From Eqs. (fT9")l and (|23|) . the final expression for the 
effective potential up to order ring can be written as 
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(28) 



where we have dropped out the explicit w-independent 
terms. The potential is real provided that to 2 + lli > 
and m 2 + II* > (this last condition is enough to make 
sure that m 2 . + U° > 0). The first condition holds if it 
happens for v = which yields the requirement that 
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FIG. 5: The second derivative of V e // at v = as a function 
of c and T for eB = 0. 
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The second condition defines a a magnetic field depen- 
dent value of a temperature Tb below which the analysis 
breaks down. 



IV. RESULTS 

To explore the properties of the effective potential 
in Eq. I|28p. we first study the case with B = and 
c = 0.25. Figure [3] shows the effective potential scaled by 
Ti(0.25) 4 , dropping out u-indcpcndent terms, computed 
for the same values of the parameters as in Eq. I|20p . No- 
tice that for the highest temperature the minimum of the 
effective potential is at v = 0. As the temperature de- 
creases, the potential flattens at v = with its curvature 
remaining positive. However an interesting phenomenon 
happens, namely, that for even lower temperatures, the 
potential becomes steeper at v = 0. This behavior per- 
sists down to the lowest temperature allowed. This effect 
is caused by the ring diagrams even in the absence of the 
magnetic field and thus precludes the phase transition to 
occur. 

To see how this behavior is affected by varying the 
model parameters, we keep B = but change the pa- 
rameter c which we now take as c = 0.38. Figure [4] shows 
the effective potential scaled by Ti(0.38) 4 , dropping out 
w-independent terms, computed for the same values of 
the parameters as in Eq. ([2U)l . Notice that for the highest 
temperature the minimum of the effective potential is at 
v = 0. As the temperature decreases, the potential flat- 
tens and the curvature at v = becomes negative thus, 
a minimum at a finite value of v develops. This behavior 
persists down to the lowest temperature allowed. Thus 
we see that for certain values of c the phase transition is 
allowed but below a certain critical value c c , the phase 
transition does not occur. This also indicates that, as the 
temperature decreases after the development of the min- 
imum at v 7^ 0, the system goes back to the symmetry 
restored phase. 

To find the critical value for the parameter c, we can 



see whether the condition for the curvature at v — to 
change sign is satisfied for real values of T. This is shown 
in Fig. [5] where we plot the second derivative of the effec- 
tive potential at v = as a function of c and T. Notice 
that this second derivative can change sign from positive 
to negative as the temperature decreases only above the 
critical value c c ~ 0.3. Notice also from Fig. [S] that for a 
given temperature, as the value of c increases the curva- 
ture changes sign from positive to negative. This signals 
that the system starts out in the broken symmetry phase 
requiring a much larger value of the temperature to re- 
store the symmetry. 

We now proceed to analyze the behavior of the effective 
potential including the effects of the magnetic field. Fig- 
ure[6]shows the relative difference (V e f f - V^ff°)/V^ff° 
computed for the same values of the parameters as in 
Eq. (2D]) and a field intensity eB = 0.9 x (m™ c ) 2 . For 
the chosen temperatures, this relative difference amounts 
only for a 5% — 10% showing that up to this order in the 
approximation the effects of the magnetic field are small. 
The difference becomes a bit more significant for the low- 
est temperatures. 



V. CONCLUSIONS 

In conclusion, we have studied the chiral phase transi- 
tion in the linear sigma model, including quarks, at finite 
temperature in the presence of weak magnetic fields. This 
has been accomplished by looking at the effective poten- 
tial up to the contribution of the ring diagrams, which in 
the weak field limit, is the first order where the effects of 
the magnetic field become evident. 

These diagrams are in principle an important ingre- 
dient in the calculation, given that the theory contains 
massless quantum modes which make it necessary to 
consider a resummation scheme. To accomplish the re- 
summation, one needs to compute the pion self-energy. 
Nonetheless, given the fact that the theory has large cou- 
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pling constants, the computation of such self-energy is 
not perturbatively reliable. In an attempt to qualita- 
tively study the structure of the theory, we have gone 
ahead and computed this self-energy perturbatively but 
have parameterized our lack of a non-perturbative calcu- 
lation introducing a factor c which we allow to vary such 
that < c < 1. In these terms, the analysis shows that 
considering the ring diagrams, there is a critical value 
for c below which the phase transition is precluded from 
happening. Above this critical value, the phase transition 
keeps being second order but as the temperature drops 
the system comes back to the symmetry restored phase. 
This means that if the phase transition was completed, 
then the bounce back does not have any effect. However 
if the transition is delayed, part of the system could get 
trapped in the symmetric phase. Finally, the effects of 
the magnetic field, in the weak field limit, only account 
for differences up to about 5% — 10% and become larger 
as the system cools down to the lowest temperatures for 
the analysis to be valid. 

In order to be able to extract more reliable conclu- 
sions, it is clear that one needs to have a better control 
on the calculation of the self-energy for which a non- 
perturbative treatment is called for, including the effects 
of the heat bath as well as those of the magnetic field. 
Nevertheless, the analysis shows in a qualitative way that 
if this calculation happens to yield values for the self- 
energy in the vicinity of the ones parameterized by c c , 
the magnetic field could retard the phase transition. This 
line of thinking is to be pursued in the future. 
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Appendix: Computation of ring diagrams 

For the sake of the argument, consider only a single 
charged scalar field. The generalization to the sigma 
model is immediate once we consider the isospin factors. 
The effective potential at ring-order is given by 



y(ring) 



2 

y(2) 



d 3 k 
(2tt) 3 



{ln[l + n A ] — n A } 



(30) 



where is the two-loop contribution to the effective 
potential. Notice that the charged scalar Matsubara 
propagator is given as in Eq. j24[) and that we can write 
n B = III + II, where IT and II are the magnetic field in- 
dependent and dependent pieces of II s . Thus, for small 
magnetic fields, keeping only the leading terms in the 
first equation, we get 



y(ring) 



T 



(27T) 
III 



dkk 2 I In 



iF 



{eBf 



(k 2 + m 2 ) 3 



2{eB) 2 k\ 
{k 2 + m 2 ) 4 



k 2 
IF 

y( 2 : 



(eB) 2 2{eB) 2 k\ 
(fc 2 + m 2 ) 3 k 2 +m 2 ) 4 



(31) 



The logarithmic term in the above equation can be writ- 
ten as 



In 



IF 



k 2 



IF 



(eB) 2 2{eB) 2 k]_ 
{k 2 + m 2 ) 3 ~ (fc 2 + m 2 ) 4 



In 
In 
In 



1 + 
1 - 
1 



IF 



(eB) 2 U 



2k 2 



k 2 +m 2 + U B V (fc 2 + m 2 ) 2 {k 2 + m 2 ) 3 



k 2 + r 

k 2 + m 2 + IF 3 V (k 2 + m 2 ) 2 (k 2 + to 2 ) 3 



2fc 2 



(32) 



With this expansion the ring contribution to the effective 
potential looks like 
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y(ring) 



T 



(2^) 2 Jo 



dkk 2 { In 



1 



IF 



k 2 + m 2 



IF 



k 2 + to 2 



2fc 2 



1 



(eB) 2 nx 



1 



2k\ 



1 



(k 2 + m 2 ) 2 (k 2 + m 2 ) 3 J k 2 + m 2 + H B \(k 2 + m 2 ) 2 (k 2 + m 2 ) 3 J k 2 + m 2 

I 



F (2) .(33) 



In order to render an analytical manageable expression 
let us make the approximation in the above equation such 



that 



J 



2k 2 



(k 2 +m 2 ) 2 (k 2 + m 2 ) 3 



(k 2 



r 



if 8 ) 2 (k 2 



2k\ 



IF 3 



(34) 



The error involved is of the order of the leading term in 
II, that is III which could be large if this last is propor- 
tional to the coupling A, as is the case in this work. In 
theories with a small coupling this error is small. In the 
absence of a non-perturbative calculation of fix we will 
restrict ourselves to this approximation to keep track of 
the analytical structure of the effective potential. 
Thus, under this approximation, we get 



y(ring) _ J_ 



T 

4tt 



m 

T 



(m 2 + n B ) 3 / 2 mlF 



(eB) 2 ^ 
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m 3 (m 2 + n s ) 3 / 2 



V< 2 



(35) 



We now recall that the two-loop contribution to the ef- 
fective potential contains a term that can be written as 



X 



(eB) 2 T 
967rm 3 

( e g) 2 rn! 

247TTO 3 A 

m 2 T 2 



2Ui 
~J~ 

All 2 
A 2 

(24p + ~M^ 2 



%T 3 



96tt 



(eBfrn.1 

1927rm 3 



where in the last line we have used that for this theory, 
up to one-loop order we have 



rii = a 



11 

24 



(37) 



Finally, recall that the one-loop contribution to the ef- 
fective potential contains a term that looks like 



m 3 T 
' 12tt 



(38) 



By adding +V ( - rmg '> we see that the potentially dan- 
gerous terms with odd powers of m all cancel and that 
this cancellation effectively amounts for the replacement 
ni — > y/m 2 + H B . The same discussion applies in the 
(36) case of the self-energy as can easily be checked. 



[1] For a recent and comprehensive review, see the proceed- 
ings of the workshop on high pt physics at LHC (LHC07), 
Eur. Phys. J. 155 (2008). See also the proceedings of 
Quark Matter 2009, to appear in Nucl. Phys. A. 

[2] D. E. Kharzeev, L. D. McLerran and H. J. Warringa, 
Nucl. Phys. A803, 227 (2008). 

[3] Z.-T. Liang and X.-N Wang, Phys. Rev. Lett. 94, 102301 
(2005). 



[4] F. Becattini, F. Piccinini and J. Rizzo, Phys. Rev. C 77, 
024906 (2008). 

[5] V. P. Gusynin, V.A. Miransky and I. A. Shovkovy, Phys. 
Lett. B 349, 477 (1995); Phys. Rev. D 52, 4747 (1995); 
Nucl. Phys. B 462, 249 (1996); Nucl. Phys. B 563, 361 
(1999); UKr. J. Phys. 45, 603 (2000). D.-S. Lee, C. N. 
Leung and Y. J. Ng, Phys. Rev. D 55, 6504 (1997). D. 
K. Hong, Phys. Rev. D 57, 3759 (1998), E. J. Ferrer 



10 



and V. de la Incera, Phys. Lett. B 481, 287 (2000). A. 
Ayala, A. Bashir, A. Raya and E. Rojas, Phys. Rev. D 
73, 105009 (2006); E. Rojas, A. Ayala, A. Bashir and A. 
Raya, Phys. Rev. D 77, 093004 (2008) ;K. G. Klimenko 
and V. Ch. Zhukovsky, Phys. Lett. B665, 352 (2008). 

[6] A. Sanchez, A. Ayala and G. Piccinelli, Phys. Rev. D 75, 
043004 (2007); A. Ayala, A. Sanchez, G. Piccinelli and 
S. Sahu, Phys. Rev. D 71, 023004 (2005). 

[7] E. J. Ferrer and V. de la Incera, Phys. Rev. Lett. 102, 
050402 (2009). 

[8] E. S. Fraga and A. J. Mizher, Phys. Rev. D 78, 025016 
(2008). 

[9] M. E. Carrington, Phys. Rev. D 45, 2933 (1992). 
[10] N. Petropoulos, Linear sigma model at finite temperature, 



|arXiv:hep-ph/ 0402136 
[11] N. Sadooghi and K. Sohrabi Anaraki, Phys. Rev. D. 78, 
125019 (2008). 

[12] K. Fukushima, D. E. Kharzeev and H. J. Warringa, Phys. 

Rev. D 78, 074033 (2008). 
[13] N.O. Agasian and S.M. Fedorov, Phys. Lett. B663, 445 

(2008). 

[14] B. C. Tiburzi, Nucl. Phys. A 814, 74 (2008). 
[15] A. Das, Finite Temperature Field Theory (World Scien- 
tific 1997). 

[16] M. Le Bellac, Thermal Field Theory (Cambridge Univer- 
sity Press, 1996) 



